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Abstract 

We consider the chiral condensate in the baryonic dense medium using the generalized 
Sakai-Sugimoto model. It is defined as the vacuum expectation value of open Wilson line 
that is proposed to be calculated by use of the area of world-sheet instanton. We evaluate 
it in confined as well as deconfined phase. In both phases, the chiral condensate has a 
minimum as a function of baryon density. In the deconfined phase, taking into account 
the chiral symmetry restoration, we classify the behavior of chiral condensate into three 
types. One can set the parameter of the theory such that the results, in low but sufficiently 
higher density, is in agreement with the expectation from QCD. 
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1. Introduction 

One of the most important problem in strong interaction physics is the understanding 
of the chiral symmetry and calculation of its order parameter in dense medium, since it is 
the key quantity determining the hadron property in the nuclei as well as in neutron stars. 
For long time, there has been many speculation on the behaviour of chiral condensation 
in dense medium: one of the intuitively compelling one is that as an order parameter, it is 
non-zero at vacuum and vanishes at the large enough density where the chiral symmetry 
is restored. Therefore the most naive but natural behaviour is to decrease as a function 
of density until it vanishes at the transition point [1]. However, to our knowledge, there 
is no work which prove this scenario from the first principle. This is because the strongly 
interacting nature and the presence of the chemical potential have been blocking any 
reliable calculation even in the numerical approach: the Dyson- Schwinger equation for the 
resummation is not justified in strong coupling since it involve the truncation of significant 
part of the diagrams. Furthermore the solution of notorious sign problem in lattice is not 
available yet. 

Since this situation, there have been many activities to utilize the gauge/gravity cor- 
respondence [2] to improve the understanding of the nuclear physics. Most notable one is 
the model of Sakai-Sugimoto [3,4], since the chiral symmetry and its breaking is intuitively 
realized in geometric way. However, it has its own drawback, since one can not include 
the quark mass in natural way and calculating chiral condensation is blocked for similar 
reason. So far, in order to improve this problem, many prescriptions (e.g. the uses of 
tachyonic DBI action [5-8], the introduction of additional D4-branes or D6-branes [9,10], 
etc.) have been suggested. Especially in this paper we are interested in the idea suggested 
by Aharony and Kutasov [11] to utilize the world-sheet instanton for quark mass and chiral 
condensation. The idea is to separate the ip and ip such that one is in D8-branes and the 
other is in anti-D8-branes in the UV region and covariantize the composite operator by 
inserting the open Wilson line connecting them. Then the gravity dual of this operator is 
the world-sheet instanton, a minimal area whose boundary is given by the open Wilson 
line connecting i/S and ip. This is analogous to the Wilson loop calculation in the closed 
Wilson line case [12]. 

The baryons in the Sakai-Sugimoto model are realised as the instantons, in other 
words, Skyrme-like solitons from the viewpoint of the flavour D8-branes [13], while the 
baryons are also interpreted to the baryon vertex given by the D4-branes wrapping S 4 , 
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which we call baryonic D4-branes [14]. Therefore one can describe the baryonic medium 
by the use of the baryonic D4-branes which are uniformly distributed in M 3 [15]. The 
Sakai-Sugimoto model in dense medium has been attracting many interests [15-18] as well 
as the other holographic models have [19,20]. In a recent paper, the authors of Ref. [21] 
worked out the problem of chiral symmetry order parameter in baryonic matter in the 
context of holographic NJL model [22], which is uncompatified version of Sakai-Sugimoto 
model, where there is no confinement. In this paper we will work out the calculation in 
compactified version with confinement. 

We will find that we can tune the parameters of the Sakai-Sugimoto model such that 
chiral condensation decreases as function of density until density reach the chiral transition 
point, where it suddenly drops to zero. That is, it does not decreases to zero continuously 
but will go to the first order transition. 

The remaining of the paper goes as follows: In Section 2, we review the definition 
of the generalised Sakai-Sugimoto model [23] and the open Wilson line operator as chiral 
condensate [11]. In Section 3, we analyse the confined phase. After showing the force 
balance condition by following Ref. [15], we numerically calculate the chiral condensate, 
and compare it with the normal nuclear density by the use of the experimental values. In 
Section 4, in the way similar to the confined phase, we investigate the chiral condensate 
and the chiral symmetry restoration in the deconfmed phase. Section 5 is devoted for the 
conclusion and discussion. 

2. Chiral condensate in generalised Sakai-Sugimoto model 

2.1. Generalised Sakai-Sugimoto model 

The Sakai-Sugimoto model [3,4] consists of N c colour D4-branes and Nf pairs of flavour 
D8-branes and anti-D8-branes. The open strings connecting the D4-branes and the D8- 
branes (anti-D8-branes) describe quarks (anti-quarks). In large N c where the AdS/CFT 
correspondence is valid, we can interpret the colour D4-branes as a background which has 
two phases. 
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fig. 1 (a) The confined phase where the chiral symmetry is broken. 

(b) The deconfined phase where the chiral symmetry is broken. 

(c) The deconfined phase where the chiral symmetry is restored. 



One is the confined phase (fig. la), the other is the deconfined phase (fig. lb,c). The 
transition between these two phases is realised as Hawking-Page transition. We can deal 
with the flavour D8-branes and anti-D8-branes as a probe in this background. 
The background (Euclidean) metric in the confined phase is described as 



R~ 2 ds 2 
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Note that the coordinates are dimensionless due to rescaling by the radius R which has a 
dimension of length. O4 is the volume of a unit S 4 , i.e., O4 = 87r 2 /3. The D4-branes are 
originally spanned in the (t, xi, x%, X3, X4) directions where t is Wick-rotated as t — > ir. 
The X4 direction is compact with a period ^4, namely, x^ ~ X4 + Although generally 
this yields a conical singularity at u = «kk, one can remove it by choosing 
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so that the u-x<± space has cigar geometry (see fig. la). The period (3^ leads to the Kaluza- 
Klein Mass and the Yang-Mills coupling, 
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Therefore one can obtain the relations between the parameters in the gravity side, R, «kk 
and g S: and those in the Yang-Mills theory side, Mkk and A (:= gy M N c ), as follows: 

In order to introduce a temperature in this phase, we compactify the r direction with a 
period /3 T , so that the temperature is given by 2n/(R/3 T ). 

In the deconfined phase (fig. lb,c), the background metric is 



R~ 2 ds 2 = u* 

i=l 



3 



f T (u)dr 2 + J2(dxi) 2 + (dx 4 f 



_ 3 
+ U 2 



dv? 



+ u 2 dVt\ 



(2.4a) 



e* = g s ui, F 4 = (27r) Y )liVc e 4 , R 3 = ng s N c (a')l , (2.46) 

/t(«) = 1-(^) • (2.4c) 

This Euclidean metric allows us to consider thermodynamics in the deconfined phase. The 
r direction is compactified with a period (3 T , i.e., r ~ r + (3 T , so that the period should be 
related to ut as 

by the same reason as (2.2). Then the Hawking temperature T is denoted by T = 
2ir/(R/3 T ). By this compactification the r-u space of the background becomes a cigar 
(fig. lb,c). We also compactify the X4 direction with the period ^4. 

In the backgrounds, (2.1) and (2.4), we consider the flavour D8-branes and anti-D8- 
branes as a probe. In the confined phase, the D8-branes and anti-D8-branes are connected 
and become U-shape D8-branes. Since the pair of D8-branes and anti-D8-branes describe 
U(Nf)i, x U(Nf)n chiral symmetry, the U-shape D8-branes imply the chiral symmetry 
breaking, U(Nf)i l x U(Nf)n — > U(Nf). On the other hand, there are two possibilities in 
the deconfined phase. One is the chiral symmetry broken phase which is represented by 
the U-shape D8-branes (fig. lb). The other is the chiral symmetry restored phase given 
by the configuration of parallel D8-branes and anti-D8-branes (fig. lc). 

In the original Sakai-Sugimoto model, the D8-branes and anti-D8-branes are placed at 
the antipodal position in the circle of X4 direction, in other words, the separation i between 
the D8-branes at the UV limit, u = 00, is a half of the period of x 4 direction, £ = /3 4 /2. 
However one can generalise the separation to be any £ less than ($4/2. In this paper we 
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consider this generalised Sakai-Sugimoto model, to which Bergman et al. introduced a 
baryon density [15]. 
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fig. 2 The V-shape D8-branes with baryonic D4-branes. The green shaded region 
is the world-sheet of string corresponding to the open Wilson line. 

The baryons are realised as the D4-branes wrapping S which are located at the tip 
of the flavour D8-branes. Therefore the system of generalised Sakai-Sugimoto model with 
baryons should obey the balance condition for the forces between the D4-branes and the 
D8-branes. This condition gives rise to the modification of the D8-branes from U-shape 
to V-shape (fig. 2) [15]. 

2.2. Open Wilson line as chiral condensate 

So far many trials for incorporating a chiral condensate into the Sakai-Sugimoto model 
have been done. In this paper we adopt the method introduced by Aharony and Kutasov 
[11]. Let us consider the open Wilson line operator, 
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where $ is one of the scalar fields in the super- Yang-Mills theory and i,j are indices 
of the fundamental representation of U(Nf). This corresponds to the well-known order 
parameter of chiral symmetry breaking, ip^t/jRi, in QCD. In the D-brane configuration 
that we are considering, the operator ip is at the UV boundary on the V-shape (or U- 
shape) D8-branes. Therefore, following the AdS/CFT correspondence, one can calculate 
the vacuum expectation value of 0\ , 
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where Sq is the on-shell Euclidean action of the open string whose world-sheet is bounded 
by the flavour D8-branes (see fig. 2). To leading order in a/, Sq is calculated by the 
minimal area of the string world-sheet, 

S = ^ ^d 2 ay/feUj, (2.6) 

where S is the region surrounded by the flavour D8-branes (see fig. 2). (O) is proportional 
to m q {ipip) , and indeed Ref. [11] has shown that (O) naively satisfies the Gell-Mann-Oakes- 
Renner relation, vn\j 2 oc (O). 



3. Confined phase 

3.1. D-brane configuration and force balance condition 

In this subsection, following Ref. [15], we shall show the V-shape solution of flavour 
D8-branes, the baryonic D4-branes wrapping S 4 and their force balance condition. 

We start with the action of the Nf flavour D8-branes embedded in the background 
(2.1) which have the world-volume coordinates (t,xi,X2,X3,u,Q,4) and the collective co- 
ordinate x^(u). The world- volume gauge fields on the D8-branes are decomposed into 

A = A su(Nf) + -j== A u{1) . 

Since the £7(1) part is related with the baryon density that we are here interested in, we 
turn on 

a Au) ■= rjp^=A u(1)jT (u) , (3.1) 

where we assumed that the gauge field depends only on u. Note that the rescaled gauge 
field a T (u) is dimensionless. Then the DBI action of D8-branes is written down as 

S (DBI) = _^y s r drduC ^ ^ ( 3 _ 2a ) 



CWtM = n^f(u){x' 4 (u)) 2 - ^ (a») 2 + , (3.26) 

where Af = Nf/jgQiR 9 /g s and V^(= J d s x) is the volume of M 3 . fj>$ is the D8-brane's 
tension, i.e., // 8 = (2n)~ 8 (a')~ 9 / 2 . Note that V 3 is a dimensionless value because the 



coordinates of M 3 , Xi (i = 1, 2, 3), are dimensionless. We introduce an electric displacement 
field d(u) by 



d{u) := 
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= ua' T (u) 



f(umu)) 2 -±(aUu)) 2 + ^ 



(3.3) 



5a' T {u) 

This will later be associated with the baryon density. In terms of d(u), the equations of 
motion for x±(u) and a T (u) from the action (3.2) are expressed as 
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du 
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Integrating these equations once with respect to u, we obtain 



(x',(u)) 2 = 



1 



f(u)(u s + u 3 d 2 ) 

/K)K + ^ 2 ) 
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d(it) = (i , 



(3.4a) 
(3.46) 

(3.5a) 
(3.56) 



where d and uq are integration constants, uq is the point at which x' 4 diverges (see fig. 2). 
The Chern-Simons action of the D8-branes, 

Sf s > = f /ft A TV(2WJ") 3 = ^ / c 5 (^) . 
where u 5 (A) = T^AF 2 — 2~ 1 „4 3 J r +10 _1 .4 5 ), introduces a source of a T (u) at u = u c (> w ): 



A^ C F 3 / drdu ■ 
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We assume a uniform distribution of baryons in the R 3 of x l (i = 1, 2, 3), namely, 



1 



rtri^ 



(3.6) 



ri4 implies a baryon density and indeed one can measure the realistic baryon number 
density ns with (length) -3 dimension by 

ri-4 



n B '■= 



R 3 



(3.7) 



TI4 also corresponds to a dimensionless density parameter of baryonic D4-branes wrapping 
S 4 . By (3.1), the source term can be written as 



Y (source) 

N C R 2 V 3 
' 8 2na' 



drdua T (u)n4d(u — u c ) . 



Therefore the equation of motion for a T (u) coming from the total action, S]? m + g^ bOUTCe ^ 
is 

N R 2 

Afd'(u) = — -n±8(u — u c ) . 
2iva 

By integrating this with respect to u, the densities, d and n^, are related by 

N C R 2 . . 

d=~ 7T7™4 • (3.8 

The source of baryon charge can be interpreted as -/V 4 (= n 4 V 3 ) baryonic D4-branes 
wrapping S" 4 at u = u c . Since the DBI action of the D4-branes is 



£4 = — N4fi 4 J drdQ.^e ^ a/ det g , 



where ^4 = (27r)- 4 (c/)- 5 / 2 , the free energy £4 of the on-shell D4-branes wrapping S 4 at 
u = u c is given by S 4 | on _ she ii =: - fdr£ 4 , 

AfV 3 

£4 = —^—u c d. (3.9) 

We now consider the force balance condition between the baryonic D4-branes and the 
flavour D8-branes. The force of the D4-branes along the u direction at u = u c is evaluated 

as 



u = 1 



d£ A AfV 3 , 4 



Since the Legendre transformation of the D8-branes' action (3.2) with respect to a T {u) 
becomes 

S 8 = 4 DBI) -AfV 3 J drdud(u)a' T (u) 

= ~ MV ° I dT />»V(«w) 2 ^)(i + ^ 

we can read the tension of D8-branes at u = u c under the on-shell condition (3.5), 

JW„ 3 

h = ^-utyfi% + #, (3.11) 

while the angle between the u axis and the D8-branes is described as 

y/g~^du 



cos 9 = 



\/g U udu 2 +g 4 4dxl 



u=u c 



f(u )(u 8 + uld 2 ) 



Therefore the force yielded by the D8-branes along the u direction is /§ cos 9. Finally, from 
(3.10), (3.11) and (3.12), the force balance condition, f 4 = / 8 cos6>, leads to 



3 V f(u c ) V f{u c ){ul + uld*))- {6A6) 

One can show that u c is related to £, the separation between the D8-branes in the x 4 
direction, by 
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£ = 2 dux' 4 (u) = 2 du^- 

Ju c Ju c U2 J 



u 2 f{u) 



f(u)(u 8 + u s d " , 



(3.14) 



Note that £ < /3 4 /2. The £ = /3 4 /2 case corresponds to the original Sakai-Sugimoto model. 

3.2. Chiral condensate 

In order to evaluate the chiral condensate, we need to compute the area of string 
instanton world-sheet (2.6), 

r 2 r 1 / 2 f°° r 2 f z i 2 f°° 

So = 7 / dx 4 I du = — - / dx 4 / du . 

2%a J_ e/2 J u ( Xi ) Tva J J u {x A ) 

Since there is UV divergence in So, regularisation is necessary. In Ref. [21] a UV cutoff 
parameter was introduced, however this regularisation neglects the contribution from 
the shape of D8-branes in the region of u > «oo. Here we regularise So by subtracting a 
density independent but infinite constant, 

R 2 f i/2 f°° 
Sao := — - / dx 4 / du , 

™ Jo Ju KK 

which is the area supported by two lines, x 4 = ±£/2. Then the regularized area Sq^ is 
given by 

R 2 R 2 f°° 

S T J; g := So - Soo = - — 7 £u K k : / duux' 4 (u) . (3.15) 

2na' iva' J Uc 

Then the vacuum expectation value of open Wilson line operator is defined by 

(O) reg :=e- S o\ (3.16) 



and it is identified as the chiral condensate (ipip). 

9 



We shall compute it numerically, since computing Sq S analytically is not easy. If we 
rescale the variables as 
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z :-- 
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then one can exclude the appearance of «kk from the numerical analysis. (3.13) leads to 

(3.18) 

and (3.14) can be rewritten as 

(z 3 - l)(z 5 + p 2 ) 



Z (z c ,p) := (z$ - l)(z 5 +P 2 ) = (4 - +P 2 )( 1 - \p 2 \^ 
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£ = 2 dz 
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Z (z c ,p) 
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(3.19) 



£(z c , p) is the monotonically decreasing function with respect to z c for any fixed p, and 
especially £(1, p) is equal to y/uKK^/2 = 27r/3, which means the antipodal position. The 
regularized area (3.15) is 
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Gconf 



n B 





p 0.16 fm - ^ 

(a) _(b) 
fig. 3 (a) The numerical plot of — G C onf (p) with £ = 1/2 fixed, 
(b) The critical density ns, ci(£)- 

From now on we set £ = 1/2 for definiteness of our model. The numerical plot of 
-G coni (p), the regularized area, is depicted by fig. 3a. From this plot, we can find that 
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(C)reg becomes minimum at p CT ps 2.058. Setting iV c = 3 and Nf = 2, we evaluate the 
baryon number density at p = p CT in terms of (3.7) and (3.8), namely, 

ns " 1458*4 ^ (3 - 21) 
Here we are interested in the low density region where the chiral condensate (~ (0) reg ) is 
expected to decreases. We fix Mkk and A by the use of the generalized Sakai-Sugimoto 
model at zero density [24,25]. Using the experimental values of the mass of p meson, 
m p w 776 MeV, and the pion decay constant, ps 93 MeV, one can obtain 

M KK ~ 496MeV, Aw 61.7. 

Note that these values are for the non-antipodal model with £ = 1/2, and different from 
those in the usual antipodal Sakai-Sugimoto model. Then the number density (3.21) at 
/9 cr ps 2.058 is evaluated to be 

nB,cr ~ 1-76 fm -3 . 

In terms of the normal nuclear density, no = 0.16 fm -3 , we obtain the ratio, nB iCr /no ~ 11, 
therefore ns,cr is sufficiently dense. 

For any other £, the behaviour of {0) reg is schematically same as the case of I = 1/2, 
that is to say, the regularized area has minimum at p CT as a function of density. Once 
we choose a model by fixing £, then p CT is determined. ub,ct calculated from p CT through 
(3.21) is the function of £, which is depicted by fig. 3b. ub,ct becomes equal to the normal 
baryon density at £ ps 1.022. 



4. Deconfined phase 

In this section we concentrate on the deconfined phase of the generalised Sakai- 
Sugimoto model with a baryon density. This phase is similar to the holographic dual 
of NJL model [22]. Since the fifth direction (x^) is not compact in the holographic NJL 
model, the Kaluza-Klein scale Mkk does not appear. On the other hand, in our case the 
X4 direction is compactified, so that we can compare our results with experimental values 
by calibrating Mkk- The open Wilson line in the holographic NJL model with a baryon 
density has been studied by Ref. [21], whose analysis is quite similar to ours in the gen- 
eralised Sakai-Sugimoto model. However we shall use the different regularisation scheme 
from Ref. [21] in order to include more accurate UV behaviour, and compute physical 
values by the use of the calibrated M K k and A. 
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4-1. D-brane configuration and force balance condition 



We turn on the U(l) gauge field defined by (3.1). Then the DBI action of Nf flavour 
D8-branes in the background (2.4) corresponding to the deconfined phase is 



^ DBI ) = -A/V 3 /3 T J duJC[x' 4 ,a' T ] 



C[x 4 , a' T ] = u 4 \l f T (u) {x' 4 {u)) 2 - -g (<(w)) 2 + -g 



(4.1a) 
(4-16) 



Note that the last term in the square root of (4.16) is different from that of (3.26) by 
the absence of 1/ f factor. We recall that (3 T is the period of r direction given by (2.5), 
therefore the temperature in the deconfined phase is denoted by T = = (3/4ir)y/ur- 
In the same way as Section 3, we define the electric displacement field by d{u) = —5C/5a' T 
{of. (3.3)), 



f T {u)(x' A {u)f - ^K(ti)) 2 + ^ 



d(u) = ua' T (u) 

Then the equations of motion for x 4 {u) and a T (u) are written down as 

\ , W 2 \^/. 2 1 



(4.2) 



du 



u 4 f T (u)x' 4 (u)[l + 



f T (u){x' 4 (u)Y + 



u 3 f T {u) 



= d'{u) 



In the deconfined phase, the following constant solution is allowed: 

x 4 (u) = ±-£, d(u) = d. 
2 



(4.3a) 
(4.36) 

(4.4) 



This describes the parallel D8-branes and anti-D8-branes whose separation is £ (see fig. lc). 
Therefore (4.4) implies that the chiral symmetry is restored. 
One can also find a non-trivial solution of (4.3): 



(x' 4 (u)) 2 = 



u 3 f T (u) 
d(u) = d . 



f T (u)(v° + u 6 d 
f T (u )(u 8 + u 3 d 2 ) 



1 



(4.5a) 
(4.56) 



This solution corresponds to the V-shape (or U-shape) D8-branes (fig. lb). From (4.5a) 
we can calculate the separation between the D8-branes at u = oo, 



/>oo /-c 

= 2 dux' 4 (u) = 2 

J u c J u c 



du 



1 



12 



f T (u)(u s + u s d 2 ) 
f T (uo)(u*+u 3 d 2 ) 



(4.6) 



The source term comes from the Chern-Simons action of the D8-branes as we explained 
in Section 3. We assume the uniform distribution of baryons in M 3 whose density is denoted 
by 77,4. By integrating the equation of motion for a T {u) which is derived from the DBI action 
and additionally the source term, the constants d and 77-4 are related by (3.8) again in the 
deconfined phase. The action of N±(= 714 V3) baryonic D4-branes wrapping £ 4 at u = u c 
is calculated as 

S 4 = -^Mldu c> /^W , (4.7) 
and the force generated by these D4-branes along the u direction is evaluated as 



, AfVsM § 



I / \ 3 ' 

I I ut 



1 + 2 V a 



' i: -d«U\ \f T {u c )\. (1.8) 



2R V 3 



We now compute the force from the flavour D8-branes. From the Legendre transformed 
action, 

S 8 = -A/V 3 /3 T J duu 4 ^ (f T (u){x' 4 (u)) 2 + ^3) (i + ^P^" 
the tension of the flavour D8-branes is computed to be 

A/V3/J 



R 

Since the angle between the D8-branes and the u axis is 



f T {u c ){u° + u 6 c d z ) 

the force coming from the D8-brane along the u direction becomes f%cos6. From (4.8), 
(4.9) and (4.10), we can obtain the force balance condition, = f$cos9, that is, 



>U. 1 . , ,\ / ( x /T(MQ)(«g + «grf 2 ) 



2 I 1 - - 3 frM) = ^ + - T^M^M I ■ (4.") 

4-2. Chiral condensate 

We redefine the coordinates by 

z:= — , z Q := — , z c := — , p := — — , I := i^/u^ , (4.12) 

u T u t u t u ^ 2 
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in the similar way to (3.17). Then (4.6) leads to 

Z (z c , p) := (zl - l)(z 5 + p 2 ) = (zl - l)(z 5 c + p 2 
and one can rewrite (4.11) as 

{z 3 -l)(z 5 + p 2 



3\2 



(l + 2zi) 



36^ z*{z* - l){z* + ?) 



1 



POO 

1=2 dz 

J Zc 



(4.13) 



- 1 



(4.14) 




(a) _ (b) 

fig. 4 (a) The plot of £ as function of z c with p = 1. 
(b) The p dependence of ^ max (p). 

(4.14) includes three variables, £, p and z c . As Ref. [15] has pointed out, 1 for any combina- 
tion of £ and p, z c does not necessarily exist, because the right hand side of (4.14) with a 
fixed p is the function of z c , which has a maximum, £ max (see fig. 4a). This is a remarkable 
feature different from the confined phase. By changing p, we depict the p dependence of 
f max in fig. 4b. For I > £ max , there is no solution of (4.14), that is, the V-shape solution 
does not exist, but there is only the parallel solution (4.4). On the other hand, for the 
opposite case, as one can see from fig. 4a, there are two solutions in z c of (4.14) for the 
given value of £ (< £ max ) and p. These two solutions describe different configurations of 
V-shape D8-branes. By comparing the free energies for these two configurations, one can 
clarify that the larger z c is favoured. Therefore we hereafter use only the larger z c . 

Let us evaluate the open Wilson line under the constraints given by (4.13) and (4.14). 
The on-shell world-sheet action of open Wilson line (2.6) in the background (2.4) is 

R 2 rm 



So 



1 



»oo 

dxA I du 



1 Ref. [15] considered in the variables £ and d with fixed temperature. However we here use 
the rescaled variables (4.12), such that the temperature factor ut is excluded from the following 
numerical analyses. 
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Since this integration diverges by the contribution of UV region, we have to regularise it in 
the similar manner done in the previous section, that is, we subtract the infinite constant 
defined by 

Don := : / 0x4 / au ■ 



na' Jo J UT ^f T {u) 

Since this area is supported by the two parallel lines (4.4), is related to "chiral conden- 
sate" in the chiral symmetry restoration. We shall give more comments about this issue 
in Section 5. Then the regularised area becomes 



=S -S, 



R 2 



00 nt 
ira 



G 



deconf 



dz 



1 



-1 v/i-C" 3 / 

where we used (3.16) and (4.13). 



(z 3 -l)(z 5 + p 2 
Zo(z c ,p) 



(4.15a) 
, (4.156) 



G^deconf 



9 0.16 fm" 3 - 



Pcr (a) G>) 
fig. 5 (a) The numerical plot of — Gdeconf(p) with I = 1/2 fixed, 
(b) The numerical plot of riB )Cr (i). 



Let us set £ — 1/2 for example. The numerical plot of — Gdeconf is shown in fig. 5a. It 
implies that in the deconfined phase also (0) reg firstly decreases and then increases with 
respect to the density p. (0) reg has a minimum value at 



p C r ps 2.073 . 



(4.16) 



This is interpreted by the use of (3.21) to the dimensionful value, 



riB.cr = 1.771 fm 3 w lln 
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where we assumed %k = «t for simplicity. 2 

For arbitrary £, p cr is given by a function of £. In fig. 5b, we plot ns,cr(£) under the 
assumption of «kk = «t- The chiral condensate monotonically decreases in the region 
ns < nB,cr- A remarkable feature in fig. 5b is that ns )Cr is sufficiently larger than the 
normal nuclear density. 



4-3. Can the chiral symmetry be restored by the density effect? 

We have to clarify which solution of the D8-branes is favoured with respect to the 
baryon density, the parallel configuration (4.4) or the V-shape one (4.5). For this purpose 
we compare the free energies of these two configurations (cf. Ref. [23] in zero density). The 
energy of the V-shape is given by the flavour D8-branes and the baryon vertex of D4-branes. 
On the other hand, the energy of the parallel configuration consists of only the D8-branes, 
because the D4-branes of baryons disappear into the black hole. The discrepancy of the 
free energies between those configurations is 



5(V) , C (V) ^(parallel) 
°8 4 ' °8 



-AfV 3 P T u^A£(p; 



where 



A£(p;I)=^ dz v^ 5 + p 2 (l 



\pz c \/l 



+ - P z c 



Z c 



(z 3 -l)(zl+p 2 ) 

( 2 3_ 1)(2 5 +p2) 



dz a/ 'z 5 + p 2 



Note that there are the overall minus signs in the actions. Therefore, if AS is positive 
(negative), then the parallel configuration (the V-shape configuration) is favoured. 



2 If we consider arbitrary ukk and ut, we have to take care on the difference of scalings 
(see (3.17) and (4.12)), l( collfil , e ) = y/ ukk/u t ^(deconfme) (= £), in calibrating M K k and A. The 
decreasing behaviour of 

TlB ,cr schematically does not change. 
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p 




fig. 6 (a) The plot of A£(p; 1 = 1/2), whose zeros are p x and p x '. 

(b) The plots of p x (£) and p x >(£)- Their crossing point is approximately equal 
to (0.485,3.144). 

For example, fixing I = 1/2 again, we compute A£(p; I = 1/2). The plot for A£(p; £ = 
1/2) given in fig. 6a shows that AS has two zeros, namely, the smaller zero, p x ~ 1.571, 
and the larger one, p x > ~ 5.798. In the region, p < p x and p > p x >, the chiral symmetry 
is broken, while in the region, p x < p < p x ', the chiral symmetry is restored. In the 
I = 1/2 model, p x is smaller than p cr (see (4.16)). Therefore the chiral symmetry is 
restored before the chiral condensate (0) Teg starts to increase. Since one can regard So 
as Soo, i.e. Sp g = 0, in the parallel D8-branes solution, the chiral condensate becomes 
(0) reg = 1 in the chiral symmetry restored phase. 

Let us consider an arbitrary I. The chiral symmetry restoration/breaking points, p x 
and p x ' , do not necessarily exist for any £. The £ dependences of p x and p x > are depicted in 
fig. 6b, which says that p x and p x > do not exist when £ < 0.485. In other words, the model 
with £ < 0.485 does not have the phase of chiral symmetry restoration. On the other hand, 
in the model with I > 0.485, the chiral symmetry is restored only in the window region: 

Px < P < Px'- 




( a ) ' ( b ) 

fig. 7 (a) The plots of p cr {£), p x (£), p x >{£) and £ max (p). 

(b) The magnified plot of (a) around the crossing points of p cr , p x and p x > . 
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Combining p cr (£), p x {£)i an< ^ ^max(p) that we have calculated so far, we obtain 

fig. 7. In the region of £ > £ max , the chiral symmetry is restored, because there does 
not exist the V-shape solution but the parallel one. Also in the region of I > 0.485 and 
p x < p < p x ', the chiral symmetry is restored, because the free energy of the parallel 
solution is smaller than that of the V-shape solution. In the other region, the chiral 
symmetry is broken. From fig. 7, we can find the following three types of models on the 
chiral condensate (0) reg : 




(O) 



reg 



Per P X Px' 

(b) 



Px Per P X ' 

(c) 



-P 



fig. 8 The (regularised) chiral condensate in (a) type I, (b) type II, (c) type III. 
I. I < 0.485 

Since the chiral symmetry is always broken, in other words, the chiral symmetry is 
never restored. {0) Teg firstly decreases and then starts to increase. This is depicted 
by fig. 8a. 

II. 0.485 < I < 0.488 

In this type, p cr < p x < p x > . Therefore the chiral condensate decreases in p < p CT 
and increases in p cr < p < p x . In p x < p < p x > the chiral symmetry is restored, so 
that (O)reg is equal to one. When p is larger than p x >, the chiral symmetry is broken 
again. The behaviour of (C) r eg is depicted by fig. 8b. 

III. I > 0.488 

Since p x is smaller than p cr , before (C)reg goes up, the chiral symmetry is restored, 
namely, (0) reg becomes one. Then in p > p x > the chiral symmetry is broken again 
and (O)reg increases. This behaviour is shown in fig. 8c and this case corresponds to 
the expected result from the field theory intuition. 

In models with any £, (0) Teg almost linearly decreases in low density, and increases at very 
high density. The type II and type III have the transition to the chiral symmetry restored 
phase at p = p x . 
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5. Conclusions and discussions 



We have studied the chiral condensate given by the open Wilson line in the generalised 
Sakai-Sugimoto model with the baryon density. In order to exclude mkk and ut depen- 
dences from numerical analyses, we introduced the rescaled variables, (3.17) and (4.12). 
In both of the confined and deconfined phases, the chiral condensate firstly decreases and 
then increases with respect to the density p (equivalently We have calculated the 

critical density ng ]Cr at the turning point, which depends on £. As a result, UB )Cr is the 
decreasing function of £. For instance, in the model of £ = 1/2, n B)Cr is about eleven times 
as large as the normal nuclear density n . Therefore the range, ng < JiB iCn in which the 
chiral condensate decreases, is sufficiently large, and this decreasing behaviour agrees with 
our intuition from ordinary QCD. However the increasing behaviour in very high density 
is different from the expectation from QCD. 

There is not chiral symmetry restoration in the confined phase. On the other hand, 
in the deconfined phase, there is the region of £ and p in which the chiral symmetry is 
restored, so that we have found the three types of behaviour of {0) reg with respect to £ 
(see fig. 8). In any types, the decreasing behaviour of (0} reg in low density is consistent 
with the results which we expect from the ordinary QCD. In the models of type II and III, 
there is a transition to the chiral symmetry restoration at p = p x , especially the models in 
the type III, i.e. £ > 0.488, are in good agreement with QCD. However the chiral symmetry 
restoration occurs before (0) reg reaches one, that is to say, this transition is in first order. 

In all types the chiral condensate grows in very high density, and this behaviour 
disagrees with our intuition from QCD. A possible problem in the high density is a back- 
reaction of the baryonic D4-branes. Since the high density means the large number of 
baryonic D4-branes, the background geometry should be modified by these D4-branes. 
Ref. [26] has mentioned about such kind of modification. The most honest way to deal 
with dense baryons in the generalised Sakai-Sugimoto model is to find a classical solution 
in supergravity with N c colour D4-branes and A4 baryonic D4-branes and to deal with the 
flavour D8-branes as a probe, but this would be difficult. 

As we studied so far, the chiral condensate (O) itself defined by the open Wilson line 
is not appropriate for the order parameter of chiral symmetry. However the positivity 
of — Sq S by definition, namely, (0) reg > 1, is a good tendency to the order parameter. 
Therefore we suggest the chiral symmetry order parameter x defined by 



X := (0) veg - 1 
19 



(5.1) 



so that x vanishes only when the chiral symmetry is restored: Sq is equal to in the 
parallel D8-branes and anti-D8-branes configuration. When the chiral symmetry is broken, 
X is always positive (see also fig. 8). 
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